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INTRODUCTION 

Ionic  solvation  has  been  the  subject  of  current  interest  for  at  least  a  century.  *  Modem 
experimental  techniques^  and  the  availability  of  powerful  computers^  have  produced  real  progress 
in  the  understanding  of  this  complex  problem. 

It  is  clear  that  in  many  problems  the  knowledge  of  the  very  complex  details  of  the  close  range 
interactions  are  more  than  what  one  needs.  In  the  interpretation  of  thernxxlynamic  and  kincuc 
parameters  averages  over  large  domains  arc  predominant,  and  therefore,  simple,  yet  controlled 
theories  not  only  provide  a  way  of  fitting  data,  but  also  useful  physical  insight. 

The  success  of  the  Marcus  theory  of  electron  transfer  kinetics'*  in  solution  is  a  good  example. 
In  this  theory  the  polarization  of  the  solvent  is  evaluated  from  a  continuum  naodcl  for  the  solvent, 
that  is.  the  structure  of  the  solvent  is  completely  ignored.  The  recent  rather  spectacular  calculations 
of  Chandler  and  his  group^  have  conclusively  shown  that  within  certain  limitations,  this  is  not 
urueasonable  provided  that  the  radius  of  the  activated  complex  is  chosen  in  a  suitable  manner. 

There  are  also  other  examples  where  this  is  true.  It  was  used  by  rrtany  investigators  in  the  past,  but 
in  an  entirely  empirical  way .*5 

In  a  series  of  recent  papers^'*®  we  have  discussed  the  use  of  the  Mean  Spherical 
Approximation  (MSA)  in  combination  with  a  very  simple  asymmetric  force  model,  which  has  a 
spherical  repulsive  core.  The  excellent  correlations  for  the  thermodynamic  and  kinetic  parameters 
for  a  large  number  of  solvents  are  more  than  mere  coincidence.  We  believe  that  they  arc  due  to  the 
fact  that  the  MSA  is  not  only  a  simple  theory,  but  it  also  has  some  remarkable  variational 
properties,  that  were  discovered  quite  recently.  Liquid  state  theories  like  the  MSA  and  the  HNC 
can  be  derived  as  variational  problems  of  the  free  energy  functional,  which  is  written  in  terms  of 
the  Omstein-2iemike direct  correlation  function.  This  view  was  introduced  by  Rosenfeld**'^^  for 
systems  of  hard  objects  in  general,  and  hard  spheres  in  particular.  In  that  work  the  relation 
between  the  grand  potential,  the  direct  correlation  function  and  die  scaled  particle  theory  was 
clearly  shown.  This  analysis  was  extended  to  the  case  of  charged  hard  spheres  in  general,  and  to 
the  primitive  model  of  electrolytes  in  particular.  It  can  be  shown  that  the  thermodynamics  and 


the  pair  correlation  functions  can  be  derived  from  a  simple  rntxiel  in  which  the  ions  arc  replaced  by 
charged  shells  and  the  dipoles  by  charged  dipolar  shells.  We  have  shown  that  the  direct  correlation 
function  in  the  MSA  is  just  the  electrostatic  energy  of  the  shells.  Solving  the  MSA  is  thus 
equivalent  to  an  electrostatics  problem. 

In  the  limit  of  high  density  and  charge  the  Gibbs  energy  and  the  internal  energy  should 
diverge  to  the  same  order  in  the  coupling  parameter  (charge/iemperature  ratio)  while  the  entropy 
should  diverge  at  a  slower  rate  and  therefore  the  Gibbs  energy  and  the  internal  energy  coincide  in 
this  limit.  The  mean  spherical  approximation  (MSA)  and  the  hypemetted  chain  approximation 
(NHC)  satisfy  these  bounds,  discovered  by  Onsager  some  fifty  years  ago.  (This  is  not  true  for 
the  Debye  Hiickel  theory,  where  the  entropy  also  diverges).  Therefore,  in  this  limit,  the  HNC  and 
the  MSA  Gibbs  energy,  and  its  derivatives,  coincide.  This  is  a  very  gratifying  feature,  because  the 
HNC,  which  is  the  more  accurate  theory,  is  in  general  difficult  to  solve  numerically,  while  the 
MSA  is  analytical  in  most  cases,  and  of  a  rather  simple  form. 

In  the  asymptotic  limit  the  excess  eJectrostaiic  energy  is  identical  to  the  exact  Onsager  lower 
bound,  which  is  achieved  by  immersing  the  entire  hard  core  system  in  an  infinite  neutral  and 
perfectly  conducting  (liquid  metal)  fluid. In  the  case  of  the  primitive  model  of  electrolytes,  this  is 
equivalent  to  smearing  the  charge  on  the  surface  of  the  ion.  It  was  also  shown^^d6  (y,at  the  direct 
correlation  function  is  equal  to  the  interaction  potential  of  two  "smeared”  ions  of  radius  Rj,  Rj 
where  Rj  is  the  radius  of  ion  i  plus  some  screening  length. 

The  Onsager  process  of  introducing  the  infinite  conductor,  naturally  decouples  all  the 
different  components  in  the  system  which  may  differ  in  size,  shape,  charge  distribution  and  relative 
orientation  in  space.  As  a  result,  the  variational  free  energy  functional  in  the  high  coupling  limit 
diagonalizes,  and  the  mathematical  solution  of  the  asymptotic  problem  is  given  in  term  of  the 
geometrical  properties  of  the  individual  panicles  in  the  system. 

The  MSA  Gibbs  energy  provides  an  appropriate  functional  from  which  the  charge  pan  of  the 
direct  correlation  function  is  obtained  by  functional  differentiation. 


There  is  an  actual  numerical  test  case  that  proves  the  electrostatic  analogy  of  the  MSA  Die 
very  high  charge  and  density  limit  for  the  one  component  plasma  (poin»  ions  in  a  neutrali/ang 
background)  was  investigated  using  the  HNC  equation  by  Ng’^.  in  a  tour  dc  lourcer  nunKneal 
calculation.  The  same  results  were  obtained  by  Rosenfeld  ei  al.  using  die  analytical  solution  cf  the 
MSA.‘« 

A  second  point  is  the  use  of  spherical  cores  to  represent  non  spherical  objects.  This  also 
depends  on  the  properties  of  interest.  For  the  equation  of  state,  it  works  even  for  very  non 
spherical  objects.  There  are  numerous  examples  of  spherical! zation  in  the  literature  of  molecular 
fluids.  A  notable  example  is  given  in  the  recent  work  of  Williams  et  al.*^  The  median  diameter  of 
a  non  spherical  object  yields  a  surprisingly  good  equation  of  state,  not  only  for  hard  dumbbells  but 
also  for  real  gases  at  extremely  high  pressures.^^*^*  Other  examples  are  found  in  the  extensive 
studies  of  Nezbeda  and  his  group^^-^^  on  the  structure  of  molecular  liquids. 

The  analytical  solution  of  the  MSA  for  mixtures  of  ions  and  dipoles  with  a  non  spherical 
sticky  potential  was  obtained  some  time  ago^'^'^^  for  the  case  of  a  sphere  with  a  point  dipole  and  a 
sticky  potential  of  dipolar  symmetry,  which  represents  the  specific  non  electrostatic  interactions. 

In  this  work  it  was  shown  that  the  contribution  of  the  sticky  potential  was  merely  a  shift  in  the 
polarization  parameter,  X.  This  work  has  been  continued  to  include  short  ranged  potentials  of 
higher  (tetrahedral)  symmetry.^^ 

In  this  spirit  we  have  recently  analyzed  the  solvation  thermodynamics  of  monoatomic 
monovalent  ions  in  water  and  in  a  number  of  polar  solvents. The  basic  idea  comes  from  the 
observation  that  to  a  first  approximation  the  corrections  due  to  the  molecular  nature  of  the  solvent 
result  in  a  change  in  the  screening  length  from  1/r  to  1/r  +  XOi/r  Os,  the  shift  being  different  for 
each  ion,  according  to  its  size;  because  of  geometrical  packing  considerations  the  number  of 
solvent  molecules  that  can  be  the  first  neighbors  of  a  given  ion,  i,  vary  according  to  the  ratio  OiJOs- 
Consider  now  a  realistic  model  of  the  solvent,  with  a  soft  repulsive  core,  and  of  non-spherical 
shape.  Then,  the  distance  of  closest  approach  of  a  certain  ion  should  be  very  different  for  the 
anion-solvent  pair  and  the  cation-solvent  pair,  even  if  the  unsolvated  anion  and  cation  are  of  the 


same  size.  The  correct  statical  mechanical  model  for  such  a  situation  is  the  mixture  with  non¬ 
additive  cores  discussed  in  the  literature.^^  A  simple  way  of  accounting  for  this  is  to  assume  dial  at 
low  concentration  the  effective  radius  of  the  solvent  is 

Tsi  =  fs  +  Asi  (1) 


The  Model 

In  the  MSA  for  a  system  of  spherical  ions  and  dipoles  of  arbitrary  sizc^^,  the  Gibbs  energy 
of  solvation  in  the  limit  of  infinite  dilution  is^* 


Gid  =  -  (\  -  -"I  — 

Cs 


)  ri+5« 


(2) 


8rt  £o  V 

Here,  Zi  is  the  valence  of  the  ion,  Cq  .  the  fundamental  charge,  £$,  the  static  dielectric  constant  of 
the  solvent,  Eq,  the  permittivity  of  free  space,  rj,  the  radius  of  the  ion,  and  No.  the  Avogadro 
constant.  6®  is  equal  to  rs/X,  where  Tj  is  the  radius  of  the  solvent  (the  effective  radius)  and  the 
MSA  polarization  parameter.  X  is  calculated  from  the  dielectric  constant  of  the  pure  solvent  Cj 
using  the  Wertheim  relationship^^ 


\2(Ul)4  =  16  Es  (3) 

Although  this  value  of  the  polarization  parameter  may  be  appropriate  for  estimating  dipole-dipole 
interaction  energy,  it  gives  poor  estimates  of  the  Gibbs  solvation  energy  when  this  is  assumed  to 
equal  Gid.2»10*28  Therefore,  and  in  view  of  the  above  discussion  we  change  eq.  (2)  to 


where 


Gid(C+) 


No  (zieo)2  _  K  1 

8;tEo  V  esj  (ri+5t) 


6c 
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and  for  cations,  and  similarly,  for  monovalent  monoatomic  anions 

Gid(A-)  =  -  A  -  — — 

8k£o  V  (H+Sa) 


(4) 

(5) 


(6) 


where 


(7) 


c  Ts+Asa 

This  separation  and  paraineterization  produces  a  much  more  sensible  and  physical  model  of  Gjd  few 
ions  in  a  number  of  solvents. 

The  second  contribution  to  the  Gibbs  solution  energy  comes  from  repulsive  dipole-dipole 
interactions75.33  and  is  given  by  the  equation 


c  =  (2ieo)^(£s-l)^  (OS 

‘'"■64rteoes(ri-^5.)2  [ 

This  expression  may  be  simplified  considerably  when  one  considers  the  range  of  values  typical  for 
X.  For  water  whose  dielectric  constant  is  78.3  at  25  °C,  X  is  equal  to  2.65.  The  corresponding 
value  of  the  ratio  X(X-t-3)/2(X+l)2  which  appears  in  the  denominator  of  the  term  in  square  brackets 
in  eq.  (8)  is  0.56.  This  is  negligible  with  respect  to  the  value  of  Cs  and  may  be  neglected.  The 
ratio  appearing  in  the  numerator  of  this  term,  namely,  (3X+2)/(X+l)  is  equal  to  2.7.  Thus,  the 
expression  for  Gdd  may  be  rewritten  as 


The  simplified  expression  has  the  advantage  that  the  parameter  X  does  not  appear  on  its  own  but 
only  in  5j.  This  fact  is  important  in  applying  the  expression  to  experimental  data.  Finally,  on  the 
basis  of  previous  work^S,  only  a  fraction  0  of  the  dipole-dipole  term  is  used  in  estimating  the 
Gibbs  energy  of  solvation  so  that  the  final  expression  for  Gs  becomes 

Gs  =  Gid  +  0  Gdd  (10) 

The  expressions  for  the  corresponding  entropies  are  obtained  from  the  temperature 
derivatives  of  the  above  Gibbs  energies.  For  the  ion-dipole  term,  the  result  for  a  cation  is 
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Sid(c+)  =  —i—  -  — ^ —  A  -  M 

83t£o  [es^  (ri+5c)  (n+Sc)^  i  e^j 


1_\  d5c 
.  dT 


where  dSc/dT  is  the  temperature  derivative  of  8c  which  is  obtained  from  the  temperature  derivative 
of  Xc-  Comparing  eqs.  (4)  and  (1  i),  it  is  easily  shown  that 


SidCO) 


e-c^) 


des  U  eoG.d2(C-^)  dSt 


£^2  dT  No(ziCo)2  dT 


The  corresponding  expression  for  an  anion  is 


Sid(A-) 


e-s) 


=  Gid  (A-)  ^  8n  eoGid^(A-)  dSA 
^sj  ^  No(zieo)2  dT 


From  these  equations  it  follows  that  two  additional  parameters  are  required  to  estimate  the  entropy 
of  solvation  of  a  salt,  namely,  the  temperature  derivatives  of  the  polarization  parameters  for  the 
cation  and  anion. '0.28 

From  the  expression  for  Gdd.  one  obtains  the  following  equation  for  the  corresponding 
entropy: 


Sdd  =  Gdd[^  /  ^1  -  ff/  (“Iri  +  2.7  8i)  -  (n  5,)j  (14) 

The  derivative  dX/dT  may  be  found  from  the  Wenheim  equation  from  which  one  obtains 

dX  rX(l+X)1  1  des 

“T  '  [  3X+1  J  2£,  dT  05) 

Assuming  that  the  fraction  0  is  independent  of  temperature,  the  resulting  expression  for  the 
entropy  of  solvation  is 

Ss=Sid  +  eSdd  (16) 

This  TiKxiel  is  now  examined  with  respect  to  data  for  monoatomic  monovalent  ion  solvation  in 


water. 


Results  and  Discussion 


In  order  to  apply  the  above  model  to  experimental  data  one  must  choose  radii  for  the  ions  and 
an  extrathermodynamic  assumption  for  extracting  single  ion  solvation  energies  from  experimental 
values  for  the  salts,  namely,  the  alkali  metal  halides.  Two  sets  of  ionic  radii  give  good  fits  between 
experiment  and  theory namely,  the  values  given  by  Pauling^  and  those  extracted  from  neutron 
and  X-ray  diffraction  experiments.^^  Pauling  values  are  used  here  simply  because  they  are 
more  familiar.  The  extrathermodynamic  assumptions  used  to  separate  experimental  quantities  into 
those  for  cations  and  anions  were  examined  in  detail  for  aqueous  systems  by  Conway.^  Values 
of  the  Gibbs  energy  and  entropy  of  solvation  for  the  alkali  metal  cations  and  halide  anions 
estimated  by  Conway  are  summarized  in  Table  1  together  with  the  Pauling  radii  for  these  ions. 

In  our  previous  analysis  of  the  Gibbs  solvation  energy^**®,  the  contribution  of  dipole-dipole 
interactions  was  ignored  ( 6  =  0  in  eq.  ( 10)).  Analysis  of  the  data  for  the  alkali  metal  halides 
showed  that  the  best  values  of  5c  and  6^  (eqs.  (4)  and  (6))  assuming  Pauling  radii  for  the  ions 
were  82.4  and  18.0  pm,  respectively.  The  quality  of  the  fit  between  theory  and  experiment  is 
really  excellent,  the  standard  deviation  for  the  calculated  values  being  1.5  kJ  mol'*,  that  is, 
significantly  less  than  1  percent.  This  result  is  illustrated  in  Fig.  1  in  terms  of  a  plot  of  -Gs  against 
l/(ri-H6i).  It  is  apparent  that  the  values  of  Gs  fit  very  well  on  a  straight  line  in  the  region  where  data 
exist  but  that  the  slope  of  the  best  straight  line  through  these  data  is  slightly  higher  (75.2  kJ  nm 
mol'*)  than  that  predicted  by  theory  (69.45  kJ  nm  mol'*).  As  a  result,  an  ion  with  infinite  radius  is 
predicted  to  have  a  Gibbs  solvation  energy  of  35.7  kJ  mol'*  when  this  simple  model  is  applied. 

This  can  be  attributed  to  the  failure  to  account  for  dipole-dipole  interactions  which  constitute  a 
positive  contribution  to  Gs. 

The  same  data  were  re-examined  on  the  basis  of  eq,  (10)  using  the  expression  for  Gdd 
given  by  eq.  (9).  The  fit  was  carried  out  so  that  the  resulting  estimates  of  Gid  and  Gdd  would  go  to 
zero  for  an  ion  of  infinite  radius.  The  resulting  best  values  of  6t  and  6a  were  much  smaller  than 
those  obtained  by  the  simpler  analysis,  and  are  recorded  in  Table  2.  The  corresponding  estimates 
of  Gid  are  plotted  against  lAq+Si)  in  Figure  2.  As  expected  these  estimates  are  larger  in  magnitude 
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than  the  experimental  values  of  Gg  also  plotted  on  the  san^c  graph.  Both  plots  are  linear  in  l/(ri+5i) 
and  have  zero  intercepts  within  the  experimental  standard  deviation.  Estimates  of  0G<jd  for  the 
same  nine  ions  are  also  plotted  in  this  figure.  For  the  cations,  this  quantity  varies  from  a  high  of 
41.2  kJ  mol'i  for  Li+  to  a  low  of  22.7  kJ  mol'*  for  Cs+,  and  is  approximately  eight  percent  of  the 
magnitude  of  Gid-  In  the  case  of  the  anions  it  represents  a  slightly  higher  fraction  of  the  total  Gibbs 
energy  and  varies  from  a  high  of  42.7  kJ  mol"'  for  F"  to  a  low  of  30.8  kJ  mol"'  for  I"  ion.  The 
values  of  6Gdd  are  also  approximately  linear  in  l/fn+S)  with  an  intercept  coiresponding  to  zero  for 
an  infinitely  large  ion.  The  value  of  5  used  in  estimating  Gdd  is  that  estimated  on  the  basis  of  the 
Wertheim  equation,  namely,  51.6  pm.  W^en  the  contributions,  Gjd  and  0Gdd  are  added  one 
obtains  estimates  of  Gs  which  agree  with  those  obtained  from  experimental  data  to  within  a  few  kJ 
mol"'.  Values  of  the  parameters  used  in  our  calculations  including  the  fraction  0  arc  summarized  in 

Table  2. 

A  successful  model  for  ionic  solvation  must  also  be  able  to  estimate  correctly  the  entropy  of 
solvation.  In  order  to  make  this  calculation  one  must  have  estimates  of  the  temperature  derivatives 
of  the  parameters  8t,  5a  and  X.  The  latter  was  calculated  using  cq.(15)  which  gives  a  value  of 
-2.48  X  10"*  K-'  for  dX/dT  on  the  basis  of  the  dielectric  properties  of  water.  Assuming  rj,  the 
radius  of  a  water  molecule,  is  equal  to  137  pm  the  corresponding  value  of  d5/dT  is  0.048  pm  K"'. 
Values  of  8c  and  5a  were  determined  from  the  single  ion  estimates  of  Ss  given  by  Conway. In 
the  case  of  cations  on  the  basis  of  eqs.  (12)  and  (16),  one  may  write 


_  ^ 
“  dT 


X 


(17) 


where 


and 


871  eoGid^(O) 
No(zieo)2 


Y  .  Ss(e)  (1  -  ^  §  *  es 


(18) 


(19) 
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ScCC")  being  the  solvation  entropy  of  the  cation.  Similar  equations  can  be  written  for  anions  using 
eqs.  (13)  and  (16).  A  plot  of  Y  against  X  for  both  ions  is  shown  in  Fig.  3.  Tliese  data  were  fit  to 
a  straight  line  passing  through  zero  (one  parameter  fit)  as  shown,  the  resulting  slope  which  equals 
both  dSs/dT  and  d5A/dT  being  0.032  pm  K'f  It  is  important  to  note  that  the  present  analysis 
suggests  that  these  temperature  coefficients  are  equal  in  contrast  to  our  previous  conclusion  based 
on  an  analysis  which  ignored  dipole-dipole  interactions.  The  estimated  error  in  the  slope  is  not 
large  being  about  4%.  It  should  also  be  noted  that  the  magnitude  of  d5c/dT  and  d5A/dT  is  about  65 
percent  of  that  found  on  the  basis  of  the  Wertheim  equation  (eq.  (15)). 

Using  the  above  coefficients  values  of  Sid  were  estimated  for  the  alkali  metal  and  halide  ions 
and  are  plotted  in  Fig.  4  together  with  estimates  of  Sdd  and  experimental  values  of  Sg  against  the 
experimental  Gibbs  solvation  energy  Gs-  It  is  clear  that  the  entropy  of  solvation  is  approximately  a 
quadratic  function  of  the  Gibbs  energy  as  suggested  by  the  MSA  model.  The  ion-dipole 
contribution  to  the  solvation  entropy  is  clearly  the  larger  varying  from  -161  J  K  'mol-i  for  Li'*'  to 
-55  J  K'*mo]'i  for  Cs'*';  in  the  case  of  the  anions,  this  contribution  is  -145  J  K'^mol  *  for  F  and 
-65  J  K'^mol'^  for  T.  The  positive  dipole-dipole  contribution  is  generally  less  than  ten  percent  of 
that  from  ion-dipole  interactions.  The  estimate  of  Ss  agrees  with  the  experimental  value  quite  w'ell 
except  for  the  heavier  halide  ions.  Thus,  in  the  case  of  f  ion,  the  estimate  of  Sj  is  higher  in 
magnitude  by  20  kJ  mol'f  This  may  be  a  result  of  some  oversimplification  used  in  the  present 
analysis,  especially  the  assumption  that  d5c/dT  and  d5A/dT  are  equal.  Nevertheless,  the  agreement 
between  theory  and  experiment  is  really  quite  good  and  demonstrates  that  the  MSA  provides  a 
simple  and  powerful  model  for  ion  solvation  in  water. 

The  present  analysis  may  easily  be  extended  to  the  other  polar  solvents  considered  in  our 
previous  work.^'28  it  is  obvious  that  the  values  of  X  for  these  solvents  that  would  be  obtained 
using  the  present  analysis  would  be  quite  different  from  those  reported  earlier.^-^S  However, 
preliminary  calculations  indicate  that  they  follow  the  same  trend.  Thus,  the  correlation  observed 
between  Xq  and  the  Gutmann  donor  number,  for  the  solvent,  and  X^.  and  the  Dimroth- 
Reichardt  acidity  parameter  should  be  maintained.  This  provides  an  important  connection 


between  the  MSA  X  and  these  empirical  solvent  parameters  which  are  often  used  by  chenasis  sn 
assessing  solvent  propertics.^^  **^ 

In  conclusion,  the  present  study  has  shown  that  vhc  dcscripuon  of  the  vanation  m  Gibbs 
energy  and  entropy  with  ion  size  is  significantly  improved  when  dipolc-dipoic  mtcracuons  for  die 
solvent  molecules  arc  taken  into  consideration.  This  is  especially  important  when  considering  the 
entropy  of  solvation  which  is  very  poorly  estimated  by  other  models.  The  simplicity  of  the  MSA 
makes  it  an  attractive  model  for  general  use  by  physical  chemists.  Funhermonr,  the  polarization 
parameter  has  practical  significance  as  a  measure  of  the  solvents’  ability  to  act  as  a  Lewis  acid  or  a 
Lewis  base.  In  future  papers,  we  will  expand  the  description  applied  here  for  water  to  other  polar 
solvents. 
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Table  I. 


Pauling  radii  and  standard  thermodynamic  parameters  of  solvation  of  nx)noatorruc 
monovalent  ions  in  water  at  298°  K 


Ion 

Pauling  Radius 

ri,  nm 

Gibbs  Solvation  energy 

Gs/kJ  moI‘* 

Solvation  Entropy 

Ss/J  K-‘  mol-‘ 

1. 

Li+ 

0.060 

-487.5 

-141 

2. 

Na+ 

0.095 

-387.6 

-109.7 

3. 

K+ 

0.133 

-314.0 

-  74.2 

4. 

Rb+ 

0.148 

-292.7 

-  62.1 

5. 

Cs+ 

0.169 

-260.5 

-  59.1 

6. 

F* 

0.136 

-456.8 

-136.8 

7. 

a- 

0.181 

-339.7 

-  79.9 

8. 

Br 

0.195 

-325.9 

-  64.4 

9. 

I- 

0.216 

-279.5 

-  41.4 

Table  2. 


MSA  Parameters  for  the  Alkali  Metal  and  Halide  Ions  in  Water  at  25  °C 


Radius  correction  for  cations,  5c/nm 

0.0692 

Radius  correction  for  anions,  5a/nm 

0.00 11 

Dipole-dipole  fraction  for  cations,  0c 

0.31 

Dipole-dipole  fraction  for  anions,  0a 

0.26 

Temperature  derivative  of  radius  correction 

for  cations  and  anions,  d5/dT/pm  K'l 

0.032 

J 


Legends 
Figure  1. 
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Figure  3. 


Figure  4. 


for  Figures 

Plot  of  the  single  ion  Gibbs  energy  of  solvation  for  the  alkali  metal  canons  (•)  and 
halide  anions  ( ♦ )  against  the  reciprocal  of  the  ion  radius,  r,.  plus  the  MSA  parameter 
5.  The  value  of  5  for  cations  was  82.4  pm  and  for  anions,  18.0  pm.  The  solid  line 
shows  the  least  squares  fit  to  these  data,  and  the  broken  line  the  value  of  the  lon- 
dipole  contribution  to  the  Gibbs  energy  according  to  cqs.  (4)  and  (6).  The  numbers 
refer  to  individual  ions  listed  in  Table  1. 

Plot  of  the  sing.'C  ion  Gibbs  energy  of  solvation,  Gs,  the  ion-dipole  contribution,  G.j, 
and  the  dipole-dipole  contribution.  Gdd.  against  the  reciprocal  of  the  ion  radius  q  plus 
the  MSA  parameter  5  for  the  alkali  metal  cations  (•,  o)  and  halide  anions  (♦,0).  The 
value  of  5  for  cations  was  69.2  pm,  and  for  anions,  l.l  pm.  The  solid  lines  show 
the  least  squares  fit  to  each  set  of  data.  The  numbers  refer  to  individual  ions  listed  in 
Table  1. 

Plot  of  the  function,  Y.  (eq.  (17))  against  the  function  X  (cq.  (18))  for  the  alkali  metal 
cations  (•)  and  halide  anions  (♦ ).  The  solid  line  was  determined  by  a  one  parameter 
least  squares  fit  going  through  the  origin.  The  numbers  refer  to  individual  ions  listed 
in  Table  1. 

Plot  of  the  single  ion  entropy  of  solvation.  5$.  the  ion-dipole  contribution,  Sid,  ^^1 
the  dipole-dipole  contribution,  Sdd.  against  the  Gibbs  energy  of  solvation  Gs  for  the 
alkali  metal  cations  (•,  o)  and  halide  anions  ( ♦ ,  0).  The  curves  are  fitted  to  data  for 
Sid  2tnd  Sdd-  The  numbers  refer  to  individual  ions  listed  in  Table  1. 
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